Searching for cavities of various densities in the Earth's crust with a low-energy v e 

/3-beam 
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We propose searching for deep underground cavities of different densities in the Earth's crust 
using a long-baseline v e disappearance experiment, realised through a low-energy /3-beam with 
highly enhanced luminosity. We focus on four real-world cases: water-filled cavities, iron-banded 
formations, heavier mineral deposits, and regions of abnormal charge accumulation that, supposedly, 
appear prior to the occurrence of an intense earthquake. The sensitivity to identify cavities attains 
confidence levels higher than 3ct and 5cr for exposures times of 3 months and 1.5 years, respectively, 
and cavity densities below 1 g cm~ 3 or above 5 g cm~ 3 , with widths greater than 200 km. We 
reconstruct the cavity density, width, and position, assuming one of them known while keeping the 
other two free, in each of the aforementioned cases. Finally, we introduce an observable to quantify 
the presence of a cavity by changing the orientation of the £> e beam. 

PACS numbers: 14.60.Lm, 14.60.Pq, 91.35.Gf, 91.35.Pn 
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Introduction - Experiments performed during the last 
decade have confirmed that neutrinos can change flavour 
with a probability that depends on neutrino energy, dis- 



tance travelled, squarcd-mass differences Am^i = m\ 



and At 



'31 



a CP-violation phase S, and 



on the mixing angles #12, 023 and #13 that describe the 
rotation between the mass and flavour eigenstates. In the 
presence of matter, coherent forward scattering modifies 
the mass-squared differences and angles, and so the value 
of this probability is altered [l|, . 

In this letter, we have used neutrino oscillations in 
matter to discover regions of under- and over-density 
compared to the average density of the Earth's crust. 
This has already been discussed in the literature, in the 
context of petroleum-filled cavities, employing either a 
superbeam [3[ or the flux of 7 Be solar neutrinos Q , and 
of electric charge accumulation in seismic faults prior to 
earthquakes [f|, employing reactor neutrinos. In con- 
trast, we have kept the cavity density as a free param- 
eter, and studied how different values affect our ability 
to locate the cavity and to determine its density. Our 
experimental arrangement considers long-baseline (1500 
km) neutrino disappearance using a low-energy (5-150 
MeV) /3-beam, which provides a high-purity sample of v e 
and allows for intense matter effects. 

Neutrino propagation in matter - The probability am- 
plitudes for the transitions D e — > Dp can be arranged 

in a column vector <J e = {ip e e fpefj, iper) 7 " which 
evolves according to id*S> e /dx = H^ e , where x is 
the distance travelled since creation and the effective 
Hamiltonian in the flavour basis is given by H (x) = 
1 / (2E V ) C/tdiag (0, Ara^,Amy U + A(x), with E v the 
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FIG. 1: (Color online) Cavity and neutrino beam with a 
fixed (a) and different orientations (b). 



neutrino energy and U the lepton mixing matrix [6[. 
The matter effects are encoded in the matrix A (x) = 
diag (-v^Gj?^ (x) ,0,0), where N e (x) = y e p (x) Na v is 
the electron number density, with p the matter density, 
TVav Avogadro's number, and y e = 0.494 in the crust 0. 
We have assumed a normal mass hierarchy (Am^ > 0), 
fixed the values of the squared-mass differences and an- 
gles to the best-fit values of Ref. Q , and the CP phase 
to zero. Our results are obtained by solving numerically 
the evolution equation described above. 

For a qualitative understanding, we will use for the 
flavour-transition probability the well-known Q for- 
mula Pl%^ = Pl%y c cos 4 0,3 + sin 4 0i 3 (with A 
A cos 2 013), where Pp c %p c is the two-flavour slab ap- 
proximation, valid on a pieccwisc constant density pro- 
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file 0, in the three- layer case that will be pertinent to 
this work. It is described by Pf v ^ 9 = | [W] n | 2 , where 
U = ^3X^2X^1, Uk = cos 4>k—i{nk'&) sin^fc, a is the vec- 
tor of Pauli matrices, and Hk = (sin 20^, 0, — cos2#j\/ fc ), 
with 0M k the matter mixing angle in the fc-th slab. The 
frequencies 4>k are given by &m M Xk/(&E V ), where Xk is 
the width of the fc-th slab, and Am^ t is the squared- 
mass difference. The numerical calculation and this ap- 
proximation are in reasonable agreement when using for 
the latter the dominant oscillation parameters Am^ and 

#12- 

Location and shape of the cavity inside the Earth — 
We have assumed the existence of a cavity of uniform 
density p C avity located within the Earth's crust, itself of 
density given by the Preliminary Reference Earth Model 
(PREM) [H. Together, the ocean, crust, and LID layers 
of the PREM have a depth of up to 80 km and an average 
density of (p®) = 3.3 g cm~ 3 . Interesting geological fea- 
tures such as porous rock cavities, mineral deposits and 
seismic faults lie in the crust and LID layers. Therefore, 
we have supposed that no part of the cavity is below 80 
km. The cavity itself has been modelled as an ellipsoid, 
and we have studied an elliptic cross section of it, with 
major axis length wq and minor axis length 0.2u>o, as 
shown in Fig. [TJi. 

Since the density profile of the Earth is radially sym- 
metric, the position of the neutrino source (S) can be 
fixed at an arbitrary point on its surface, for a given 
source-detector baseline, L n . In order to specify the lo- 
cation of the cavity, we need to set the distance do mea- 
sured along the baseline from the source to the cavity's 
surface. 

Low-energy (3-beams.- There is currently a proposal to 
use a pure, collimated beam of low-energy D e generated 
by means of the well-understood j3 decay of boosted ex- 



otic ions [11| and detected through v e 



12 



C -> 



12 



B 



[12j . Our /3-beam setup contemplates an ion storage ring 
of total length l tot = 1885 m, with two straight sections 
of length Straight = 678 m each Inside the ring, 6 He 
ions boosted up to a Lorentz factor 7 = 25 decay through 
|Hc++ -> %U +++ + e- + v e with a half-life t 1/2 = 0.8067 
s. Ion production with an ISOLDE technique [l4| is ex- 
pected to provide a rate of ion injection of g = 2 x 10 13 
s _1 for 6 He; we have assumed a highly optimistic 5000- 
fold enhancement of this value, which has not even been 
considered in the literature. 

The neutrino flux from the (3 decay of a nucleus in 
its rest frame is given by the formula [l5[ $ c . m . (E u ) = 
bElE e ^El - m\F (±Z, E e ) 6 (E e - m e ), where b = 
In 2/ (mg/ii/2), with m e the electron mass and /ti/2 = 
806.7 the comparative half-life. The energy of the emit- 
ted electron is given by E e = Q — E v , with Q = 3.5078 
McV the Q-value of the reaction, and F (±Z, E e ) the 
Fermi function. 

Based on the formalism of Ref. [l6| , we have considered 
a cylindrical detector made of carbon, of radius R = 4v5 
m and length h = 100 m, co-axial to the straight sections 
of the storage ring, and located at a distance Lq = 1500 
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FIG. 2: (Color online) Confidence levels (C.L.s) in the 
w vs. p C avity plane assuming there is no cavity. The four 
points A-D are special cases described in the text. 



km from it. The integrated number of e + at the detector, 
after an exposure time t, is calculated as 



N 



t9 ^ nh J dE »®tot(E„)Pv e ^p e (L ,E u )o-(E u ) , 

(1) 

with n ps 6.03 x 10 23 cm 3 the density of carbon nuclei 
in the detector and a the detection cross section [l2l ]. 
Since L > hot, Straight, h, we can write $ t ot {E v ) ~ 
*iab (E v , 6> = 0) (Z s traight /hot )S/(4ttL$), where $iab is 
the flux in the laboratory frame [Hj], 9 is the angle of 
emission of the neutrino with respect to the beam axis, 
and S = ttR 2 is the detector's transverse area. 

Sensitivity to cavities - Assuming that there is no cav- 
ity along the baseline Lq, we evaluate the sensitivity to 
differentiate this situation from the hypothesis that the 
neutrino beam traverses a cavity of width w, position d, 
and density p. To do this, we define 



X 2 {w,d,p) = 



[N? av (w, d, p) - Nf 

j\j no-cav 



(2) 



with iV? av (w, d, p) the number of e + , in the i-th energy 
bin, that reach the detector in the case where the beam 
traverses the cavity, and 7V"°" cav the corresponding num- 
ber in the no-cavity case. Given that 7 = 25, the neutrino 
spectrum extends from 5 to 150 MeV, in bins of 5 MeV. 
On account of the maximum energy considered, the pro- 
duction of muons, from ^ M + 12 C — > p + + X , is inhibited, 
and thus the v e — > channel is not included in this 
work. 

Before presenting the results, it is convenient to men- 
tion that, within the slab approximation, the first and 
third slabs correspond to the crust, and the second one, 
to the cavity, so that x\ — d, x 2 — w, and X3 = Lo—d—w. 



TABLE I: Minimum value of p C avity to have a 
statistical significance of 5<r or 10a C.L. with respect to 
the no-cavity case, for a given width w and exposures of 
3 months or 1.5 years. 
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Fig. [2] shows the sensitivity in the form of isocontours of 
\ 2 = la, 3cr, 5cr, and 10cr, with the cavity located at the 
center of the baseline, so that d = (Lq — w)/2 (which is 
tantamount to (f>i = (/> 3 ), and detector exposure times of 
(a) 3 months and (b) 1.5 years. Because of this relation 
between d and w, the analysis can be reduced to two 
parameters: p and w (or d). It is observed that the sen- 
sitivity improves as |p C avity — (/°e)l increases. Also, the 
statistical significance, for the cavity hypothesis, grows 
with u>, given that the cumulative matter density differ- 
ences are greater. 

In the slab approximation, the behaviour of the oscil- 
lation probability, that dictates the shape of the isocon- 
tours, is dominated by the 4>2 frequency, which is the only 
one that depends on both p and w. In fact, for a given 
energy, and since we have observed that the combinations 
of sines and cosines of (f>i in the probability vary slowly 
with u>, each isocont our approximately corresponds to a 
constant value of w^Jc^/E 2 + 2coCipcos29i2/E l/ + c\p 2 , 
where c = and ci = — Gi?y e A^Av cos 2 6*13/^2. 

Table U shows the minimum value of p C avity needed to 
reach 5<r and 10cr separations, for small (w = 50 km), 
medium (100 km), and large (250 km) cavities, after de- 
tector exposure times of 3 months and 1.5 years. 

Determination of cavity parameters.- The next step in 
the analysis is to change the assumption of no-cavity to 
that of a real cavity. This implies replacing M^°~ cliv —> 
A r f av (w , d , Cavity) in Eq- ©• Firstly, we consider a 
known cavity position (i.e., d = do, with the hypotheti- 
cal cavity no longer centered) and, heretofore, we set 1.5 
years as exposure time, unless otherwise specified, and 
keep the other details of the study equal to those pre- 
sented in the sensitivity analysis. We will study the four 
points, motivated by centered real- world cavities, marked 



in Fig. m A (p, 



cavity 



1 g cm 



0.555, w = 250 



km), corresponding to a deep water- filled underground 
cavity, B (p C avit y = 5 g cm~ 3 ,j/ e = 0.5, w = 250 km), to 
an iron banded formation [17j, C (/0 C avity = 10 g cm~ 3 , 
y e = 0.5, wo = 100 km), to a heavier mineral deposit, 
and D (/9 C avit y = 25 g cm" 3 , y e = 0.5, w = 50 km), rep- 
resenting a zone of seismic faults with the typical charge 
accumulation that supposedly exists prior to an earth- 
quake of magnitude 7 in the Richter scale [l8[ . In Ref . Q , 
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FIG. 3: (Color online) Confidence levels in the w vs. p 
plane. Plots (a)-(d) correspond to central cavities A-D 
from Fig.O The real values (po,Wo) are marked by 
crosses and the best-fit values by stars. 



a value of p C avity ~ (p©) and a maximum value of y e m 4 
at the fault are used, while, to keep in line with the anal- 
ysis of cavities A-C, we have cquivalently taken for cavity 
D Pcavity = 25 g cm -3 and y e = 0.5. 

Notice that the closest distance from the cavities to 
the Earth surface, for A and B, is about 19 km, and, for 
cavities C and D, 34 km and 39 km, respectively. Wider 
or uncentered cavities would lie closer to the surface; for 
instance, a centered cavity with wo = 323 km would lie 
only 12 km deep, which is roughly the current maximum 
drilling depth [lj|. 

Secondly, in Fig.[3]we observe that values of p close to 
(po) (equivalent to the no-cavity case), regardless of w. 
are excluded at the same statistical significance as shown 
in Fig. [2] The shape of the isocontours is explained by 
the same argument as in said figure. The uncertainty in 
uu is large in cases A and B, since in that density region 
the contours are parallel to the asymptote at p = (p©). 
In cases C and D, the uncertainty is smaller given that 
at higher densities the contour regions tend to be paral- 
lel to the horizontal axis. It is interesting to point out 
that in a real-case scenario, where there could be some 
prior knowledge about the density, it would be possible to 
constrain w significantly. For instance, for p = 1 g cm~ 3 
(case A), w = 2401 3 q km at la uncertainties, while for 
p = 25 g cm~ 3 (case D), w = 50 ± 10 km. 

In Fig. 31 we display the case where the cavity position 
is unknown, while its density is known (i.e., p = p C avity) 
and, thus, the contour regions are described by w and 
d (with w + d < Lq). In this figure, the size of the 
uncertainty in w follows from Fig. |3l Therefore, we are 
capable of determining if with reasonable precision in 
cases C and D. On the other hand, the determination of d 
improves mildly as the density increases. For example, in 
case A, d = 600±f 7 l km and, in D, d = 725±2?5 km at la 
uncertainties. Using the slab approximation, this larger 
uncertainty can be explained by the weak dependence of 
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FIG. 4: (Color online) Confidence levels in the d vs. w 
plane. Plots (a)-(d) correspond to central cavities A-D 
from Fig. [2] The real values (wq, do) are marked by 
crosses and the best-fit values by stars. 
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FIG. 5: (Color online) In solid black lines, A vs. a, for 
the four cavities A-D, together with selected 
uncertainty regions. 



the probability, for most of the relevant energy range, on 
<j)\ and </>3, and, therefore, on d. 

Searching for a cavity - We have explored the possibil- 
ity of varying the beam orientation, defined by the angle 
a measured with respect to the tangent to Earth at S, 
as a way of finding a cavity. This is shown in Fig. [TJd, 
and could be implemented by u sing a mobile neutrino 
detector, either on land or at sea [201. At an angle a, the 
beam travels a baseline L a ; for some values, it will cross 
a portion of the cavity with position d a and width w a . 
To quantify the size of the traversed portion, we have 
defined 



A (L , wq, do, a, p C avit y ) = 



7V CaV (L a ,W a ,d a , /dcavity ) 



with iV cav and iV no ~ cav the total numbers of e + between 
5 and 150 MeV, and L a , w a , d a functions of Lo, wq, do, 
and a. For this analysis, we have set the exposure time 
at t = 3 months. 



Fig. [5] shows the A vs. a curves for the four cav- 
ities A-D. The curve for cavity A lies below A = 1 
because its density is lower than (p©), while the den- 
sities of cavities B-D are higher. The Io~a to 5cta un- 
certainty regions around the curves are included, with 
a a = Ay/(l/N cav ) (1 + A) the standard deviation of A. 
Since A cav ~ 1/£q, then a a ~ L a and, given that, by 
geometrical construction, L a grows with a, this means 
that a a also grows with a, a feature that is observed in 
Fig. [5] Similarly to previous figures, the separation from 
A = 1 grows as /9 ca vity moves away from (p©). The C.L. 
achieved at the maximum deviation point is not as high 
as those shown in Fig. [21 due to the fact the analysis of A 
does not take into account the spectral shape, but only 



N" 



(La) 



(3) 



the total e + count. 



Separation from A = 1 at the 2(Ta C.L. is achieved 
for all cavities, with A and C reaching Sua, and D, 5<ta- 
Thus, the parameter A could be used as a quick estimator 
of the presence of a cavity, while a more detailed analysis, 
similar to the one performed for Fig. El could be used to 
estimate the cavity shape, i.e., its position and width, 
and, from the latter, its volume. 



Conclusions - We studied the use of a low-energy (5- 
150 MeV) /3-beam of v e , with a baseline of 1500 km and 
a large luminosity enhancement, to find the presence of 
deep underground cavities in the Earth's crust. We have 
determined the sensitivity as a function of the cavity den- 
sity p and size w (dimension of the cavity aligned with the 
neutrino beamline) , which reaches significances in the or- 
der of 5cr (3cr) for baseline-centered cavities with densities 
lower than 1 g cm~ 3 or greater than 5 g cm , exposure 
time of 1.5 years (3 months), and w greater than 200 km, 
rendering our analysis highly competitive. We analysed 
the C.L. regions of the reconstructed parameters of four 
real- world cavities in the p vs. w plane, for a known cav- 
ity position d, and also in the w vs. d plane, when p is 
known. Finally, we have considered sweeping the Earth 
in search of a cavity using an orientable neutrino beam, 
which proves to be a useful tool to detect its presence. 
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